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This paper is dedicated to Prof. Jacques Friedel, an inspirational scientist and a great man. 
His excellence and clear vision led to significant advances in theoretical physics, which spilled into 
material science and technological applications. His fundamental theoretical work on commonplace 
materials has become classic. We can think of no better tribute to Friedel than to apply a funda- 
mental analysis in his spirit to a peculiar class of materials - auxetic materials. Auxetic materials, or 
negative-Poissons-ratio materials, are important technologically and fascinating theoretically. When 
loaded by external stresses, their internal strains are governed by correlated motion of internal struc- 
tural degrees of freedom. The modelling of such materials is mainly based on ordered structures, 
despite existence of auxetic behaviour in disordered structures and the advantage in manufacturing 
disordered structures for most applications. We describe here a first-principles expression for strains 
in disordered such materials, based on insight from a family of 'iso-auxetic' structures. These are 
structures, consisting of internal structural elements, which we name 'auxetons', whose inter-element 
forces can be computed from statics alone. Iso-auxetic structures make it possible not only to iden- 
tify the mechanisms that give rise to auxeticity, but also to write down the explicit dependence 
of the strain rate on the local structure, which is valid to all auxetic materials. It is argued that 
stresses give rise to strains via two mechanisms: auxeton rotations and auxeton expansion / con- 
traction. The former depends on the stress via a local fabric tensor, which we define explicitly for 
2D systems. The latter depends on the stress via an expansion tensor. Whether a material exhibits 
auxetic behaviour or not depends on the interplay between these two fields. This description has 
two major advantages: it applies to any auxeton-based system, however disordered, and it goes 
beyond conventional elasticity theory, providing an explicit expression for general auxetic strains 
and outlining the relevant equations. 



I. Introduction by Sam F. Edwards 

At the end of the war, which had isolated France 
from the English speaking world, several French scien- 
tists moved to UK universities, in particular to study 
solid state theory. The outstanding person in the UK at 
the time was Nevill Mott in Bristol and Jacques Friedel 
moved to Bristol to work in Mott's group. At the time, 
the field of theoretical physics was moving into the use 
of field theory to elucidate elementary particle theory, a 
direction favoured in Cambridge University and in Lon- 
don. The Bristol group, however, specialised in electronic 
studies, an area that Friedel preferred. I remember his 
papers at that time, which had a wonderful clarity and 
discussed down-to-earth type of problems. It was re- 
freshingly in stark contrast to the renormalisation the- 
ory, which was the fashion in quantum field theory at 
the time. 

Sometime later, Mott moved to Cambridge and 
Jacques returned to Paris. This reminds me my first 
conference in Paris, where I gave my first paper. It was 
nonsense, I regret to say, for it tried to separate green 
functions for the real and imaginary parts of the wave 
function. Fortunately, none of the attendants in that 
conference exposed it. 

Anyway, I recall Friedel giving wonderful lectures in 
Cambridge, where his work was held at very high esteem. 
Years later Cambridge University awarded him an Hon- 
orary Doctor of Science and I had the pleasant task of 
arranging a dinner for him. Friedel was also involved in 



setting up the European physical society, where I was ac- 
tive, and I recall him giving valuable advice on its struc- 
ture. 

A central sociological problem in theoretical physics is 
to choose the problem to work, for there are many bril- 
liant people working at the forefront of the field. Think- 
ing of Friedel's work on electronics in parts of systems, it 
occurred to me that one should be able to do statistical 
mechanics on continuous systems in contrast to particu- 
late systems. With my coauthor here, Raphael Blumen- 
feld, I have developed this idea by studying the entropy 
of particulate systems in the continuum. For example, 
in conventional thermal systems the entropy S* is a func- 
tion of pressure, volume, energy and number of particles, 
S{E, P, V, N) and one of the most useful concepts it leads 
to is the temperature T = dE/dS. We have applied these 
to granular systems where the entropy is due to config- 
urational disorder and the volume takes the role of the 
energy. Consequently, the analogue of temperature is the 
'Compactivity' Xq = dV/dS. There are other quantities 
that dictate the states of granular matter, the simplest 
being the response of stresses to the entropy, X = da/dS, 
which we called the Angoricity (note that the Angoric- 
ity is in fact a tensor). An even richer and more general 
'thermodynamics' is required when we study mixtures. 

This paper is dedicated to Jacques Friedel and, in the 
spirit of the close relations of his theoretical works with 
real materials, we can think of no better tribute to him 
than to present a fundamental theory that aims to under- 
stand the physics of a peculiar class of materials - auxetic 
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materials. 

' II. General introduction 



Auxetic materials, i.e. materials with negative Pois- 
son's ratio, expand when stretched and contract when 
compressed, in contrast to most conventional materials. 
This is due to correlated degrees of freedom in the inter- 
nal elements that theses materials are made of. These 
elements are reversibly foldable and, in effect, can be 
regarded as the basic constituents of cellular solids. In 
the following, we call these foldable elements 'auxetons'. 
Macroscopic auxetic structures can be manufactured of 
polymers [Ij or metals[lj, 12|]. They can exist on a range of 
length-scales and, in particular, can be constructed out 
of molecular building blocks pL's']. Auxetic materials are 
useful in applications requiring high shear to bulk moduli 
or compactification on impact, e.g. for energy absorbing 
materials and bullet-proof armours. 

Both natural[i-[i| and man-made (Tol. [TTj auxetic ma- 
terials have been discovered, made and studied. Much 
of the theoretical analysis, however, is carried out on or- 
dered models, such as two-dimensional inverted cell hon- 
eycombs. Although models of the auxeticity phenomenon 
in ordered structures is convenient for analysis purposes, 
the ubiquity of disordered such materials and the little 
existing understanding of deformations in the presence of 
disorder require a more general theory. Here we describe 
such a theory, based on a recent suggestion made in [l^ . 

The aims of this paper are the following. First, we 
describe a new family of disordered auxetic structures, 
called iso-auxetic (lA) structures, for which it is possible 
to identify clearly the basic strain mechanisms. Second, 
we show that elasticity theory is not necessary for the 
description of auxeticity, implying that using negative 
Poisson's ratio as a descriptor has a limited utility. Third, 
we present an explicit expression for the auxetic strain in 
terms of local expansive and rotational fields. In this 
expression, the fields are coupled to the stress through 
well-defined tensors, which we discuss. Fourth, we show 
that auxeton rotations are essential to the understanding 
of the global behaviour and that the rotational field can 
be modelled without resorting to non-symmetric stresses. 
This obviates models based on Cosserat theorv[l3|. 

The paper is structured as follows. We first intro- 
duce the new family of lA structures. These are struc- 
tures whose inter-auxeton forces can be determined from 
statics alone. This property distinguishes lA from more 
conventional auxetic structures, which we term elasto- 
auxetic (EA). Specifically, the stress field equations of 
isostaticity theory differ significantly from those of con- 
ventional elasticity in that they are based on local stress- 
structure relations, as opposed to the usual stress-strain 
relations [T^ - [T6| . We next describe an extension of a re- 
cent result for yield of granular systems to lA structures 
and write down explicitly the lA strain equation in terms 
of two local fields: an expansive and a rotational. It is 
then argued that the mechanism for auxeticity depends 



only on these two fields and is therefore independent of 
the particular way that the structure transmits stresses, 
whether isostatically or elastically. Hence, the auxetic 
expression for the strain is valid for all auxetic materials. 
This, in turn, implies that general auxeticity needs to 
be described by a theory that goes beyond elasticity. A 
particular implication of this conclusion is that negative 
Poisson's ratio in auxetic materials should be regarded 
only as a descriptor of the ratio of strains in perpendic- 
ular directions, not as a ratio of elastic moduli. We also 
argue that, although the way the form of the strain ex- 
pression is the same for all auxetic materials, the strains 
developing in lA structures differ markedly from those 
developing in EA structures under the same loading con- 
ditions. We conclude with a discussion of the results. 

II. Iso-auxetic structures 

In the following discussion, we consider planar auxetic 
materials made of 2D elementary units that connect to 
their neighbours at exactly three points. We call these el- 
ements 'auxetons'. Aiming at a theory of disordered ma- 
terials, we do not require that the auxetons be identical, 
nor that the system possess any type of symmetry, trans- 
lational or otherwise. Rather, we consider systems whose 
auxetons comprise a mixture irregular sizes, shapes and 
orientations. A wide variety of such structures can be 
constructed, some of which are illustrated in figure 1. 




FIG. 1. Examples of auxetons made of three-contact building 
blocks. Each auxeton can expand and rotate when forces are 
applied to its ends, termed 'contacts' m the text. 



We constrain our auxetons to have three 'contacts' 
with their neighbours and connect these contacts by 
imaginary straight lines into triangles (the dashed blue 
lines in figure 2). This construction results in a planar 
graph of triangles, connected at their vertices. The tri- 
angles enclose polygons, which we call in the following 
cells. 

When loaded by external forces, the auxetons trans- 
mit those to one another through 'inter-auxeton' forces. 
The contacts between auxetons may or may not be free- 
jointed. One expects the latter to be more common, in 
which case a contact can support a certain threshold of 
torque moment without yielding. This gives rise to a 
finite overall stress threshold for straining the material. 
Consider then a structure, made of N{'^ 1) auxetons, 
stressed below the yield threshold by a set of external 
forces. Below the yield threshold, the system is in me- 
chanical equilibrium and all the inter-auxeton forces and 
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FIG. 2. A section of a disordered auxetic structure, made of 
joining auxetons at their contacts. The contacts are joined by 
straight lines (blue dashed) into a triangle. These triangles 
are then used to characterise the contact network in a well- 
defined manner. 



torques are balanced. Since every auxeton has three con- 
tacts then the number of contacts is 3N/2 + 0{^/N), 
where the latter term is a boundary correction, which 
can be neglected for ^ 1. Since each contact trans- 
fers one force vectors, there are overall 3A^ force compo- 
nents. These can be determined uniquely by the three 
balance equations for every auxeton - one of torque and 
two of force components. It follows that this structure 
is statically determinate, or isostatic. Hence the name 
iso-auxetic. A familiar textbook statically determinate 
system is that of a ladder on a frictional floor leaning 
against a frictionless wall. The forces that the wall and 
the floor apply to the ladder can be determined uniquely 
from its three balance equations. It is important to note 
that, as in the case of the ladder problem, the deter- 
mination of the discrete inter-auxeton forces requires no 
knowledge whatever of the elastic properties of neither 
the auxetons nor the contacts. Since the stress field is 
nothing but a continuos representation of the large num- 
ber inter-auxeton forces, it must reflect the nature of the 
discrete solution and therefore be independent of local 
elastic moduli. It follows that elasticity theory, which 
does rely on knowledge of the elastic moduli, is inappli- 
cable for I A structures. 

For later discussion, it is useful to recall the continuum 
2D stress equations of isostaticity theory - the theory of 
stresses in isostatic structures. 



aij = aji (2) 
QijCTij = (3) 

Eqs. (H)) and ^ represent force and torque balances, re- 
spectively, with a the stress tensor and g external forces. 



Eq. ([3]) is a constitutive relation between the static stress 
and the local structure, which is characterised by a sym- 
metric fabric tensor This replaces the stress- 
strain relations in conventional elasticity and is indeed 
independent of the elastic moduli of the material. 

In most known isostatic systems these equations are 
hyperbolic, leading to solutions that 'propagate' along 
characteristic paths in the material. This means that the 
response to a localised force source in 2D is generically a 
pair of force chains. In contrast, EA materials respond to 
localised sources by 'dispersing' the stress field in all di- 
rections, subject to local stress-strain relations. The dif- 
fernce between the two types of solutions stems from the 
nature of the stress field equations - while the equations of 
elasticity theory are elliptic, eqs. (ForceBal)-(ConstEq) 
of isostaticity theory are hyperbolic. The different stress 
transmission is bound to affect macroscopic behaviour, 
as will be discussed below. 

The global auxetic behaviour is the result of local fold- 
ing and unfolding of auxetons when stressed. This local 
response is independent of whether the rest of the struc- 
ture is isostatic or not, it is only dependent on the local 
expansion and rotation of the auxetons. This leads to 
the conclusion that the strain can be written in terms 
of local expansion and rotation fields regardless of the 
isostatic or elastic nature of the material. This conclu- 
sion is significant for two reasons. One is that, in lA, the 
stress is independent of elastic moduli. The other is, that 
in lA we can write the strain explicitly in terms of the 
expansion and rotation of auxetons, which means that 
the same expression holds for EA materials. This gives 
insight into the description of auxeticity in general and 
in particular into the coupling between the local strain 
and the local structure. Additionally, this suggests that 
elastic constants need not play as major a role as in con- 
ventional materials. Another important implication is 
that the negative Poisson's ratio, which such materials 
exhibit, is only a descriptor of the ratio of perpendicular 
strains and is of little use in terms of describing bulk elas- 
tic moduli because these cannot be obtained simple ho- 
mogenisation of small-scale regions. It is also worthwhile 
to note, before we continue, that this description should 
apply not only to all auxetic materials made of foldable 
auxetons, but also to those made of rigid ones [2^. 

Before we proceed, we must comment on a much de- 
bated issue: whether or not auxeticity theory necessitates 
resorting to Cosserat stress theory [IJ], which allows for 
existence of a non-symmetric stress tensor. This is not 
a question of formalism, but rather of the underlying 
physics. A symmetric stress tensor means that residual 
torque moments vanish on the continuum length-scales. 
Differently put, it means that there are no external couple 
moments on the system that require balancing mechani- 
cally by the mechanical stress field. By letting the stress 
tensor be non-symmetric on macroscopic scales, Cosserat 
theory implies that there exist external couples that the 
stress must balance. Thus, a theory that invokes only 
symmetric stresses does not resort to such additional in- 
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put and must be preferable for modelling of large-scale 
auxetic behaviour. For this reason we prefer the above 
formulation, which includes (eq. ([2]). 

Nevertheless, it is important to point out that a sym- 
metric stress tensor can still allow existence of local ro- 
tational fields of the material upon straining. In other 
words, although the stress field must be symmetric un- 
der no external couples, the strain field may have non- 
symmetric contributions. Indeed, local such contribution 
arise from rotation of auxetons and it is at the core of 
auxetic behaviour, as will be discussed in the next sec- 
tion. 

III. Auxetic strain and field equations 

To relate the strain to the local structure one has to 
have first a quantitative description of the structure, how- 
ever disordered. Such a descriptor is the fabric tensor Qij 
of eq. ([3]) . This tensor plays a key role in modelling aux- 
etic strains, as will be seen below. Consider a disordered 
structure of auxetons, comprising an arbitrary mixture 
of elements, such as those shown in figure 1. The model 
to be described below has been discussed initially in [l^ 
and it is general in that it applies to any arbitrary struc- 
ture of the above auxetons. Specifically, the disorder can 
involve both auxetons of different sizes and of different 
shapes. Connecting the three contact points around each 
auxeton by straight lines, as described above, the plane 
is tiled into a network of triangles of different sizes and 
shapes, all interconnecting at their vertices - the contact 
points. The triangles enclose polygons, which we will call 
in the following cells. According to Euler relation [2T'], a 
system of 3> 1 such auxetons would have two auxeton 
per cell. This value has small corrections from boundary 
auxetons (due to unshared contacts), but this correction 
is of order ~ 1/VN and therefore negligible. All triangle 
edge are assigned directions, making them into vectors 
r that circulate the triangles in the clockwise direction 
(figure 3). 

Every triangle is assigned a centroid, defined as the mean 
position vector of its three vertices. Similarly, every cell is 
assigned a centroid, defined as the mean position vector 
of the contacts (triangle vertices) that surround it. In 
mechanical equilibrium, the cell polygons must be convex 
to be stable. This means that, a vector R*^^ extending 
from the centroid of triangle g to the centroid of one of 
its neighbour cells c, intersects one of the triangle edge 
vectors, which we can index r'^s (figure 3). 

The vectors R'^s and r'^s can be regarded as the di- 
agonals of a quadrilateral, called 'quadron', which plays 
a significant role in granular and cellular phvsics [22H24| . 
Each quadron is associated uniquely with a pair eg (see 
figure 4). This construction allows us to quantify the 
local structure by quantifying the shape of every cg- 
quadron tensorially as the outer product 
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FIG. 3. Characterisation of the auxeton structure in 2D. We 
make the edges of the representative triangle g into vectors, 
r"^^ , by assigning the edges a direction such that they circulate 
around the triangle in the anti- clockwise direction. The vector 
IV^^ extends from the centroid of the grain contacts to the 
centroid of an adjacent cell c. 
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FIG. 4. Quadron tessellation m 2D. The vectors r*^^ and YC^^ 
make the diagonals of the cgr-quadron. The quadrons are the 
elementary units that tessellate the system. The quadron 
shape is quantified by a local structure tensor, C'i^ = r'^^R^^ . 



The tensor appearing in the isostaticity stress equations 
is the symmetric part of C"^^, summed over the cells 
around triangle g 



i=-i 

— e 
2 



E 

L c around g 



(5) 



where e is the 7r/2 rotation matrix in the plane (the Levi- 
Civita matrix) and is the transpose of C. 

Armed with a quantitative description of the local 
structure, it is possible now to relate it to the strain. 
Suppose that the structure is in mechanical equilibrium 
under a set of external forces and these forces are in- 
creased. Eventually, the system crosses what is known as 
the yield surface and it starts deforming. As will become 
clearer below, whether the deformation is auxetic or not 
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depends on the structure of the auxetons, their configu- 
ration and the magnitude of the local stresses. The aim 
of the following is to describe the equations that govern 
the strain, given the local structure and the local stress. 
Central to the model is the observation that only aux- 
eton rotation and expansion (or shrinking, which can 
be regarded as negative expansion) can give rise to dis- 
placement. The expansion corresponds to pure folding 
/ unfolding of auxetons. Thus, the local strain e, due 
to changes in shape and volumes of the triangles, can 
be written as a superposition of a rigid triangle rota- 
tion, e"""*, and triangle (non- uniform) expansion, e^^. For 
example, auxetic materials composed of rigid auxetons, 
such as those studied in can be described by e''°* 
alone. In the following we consider only the symmetrised 
strain, but there is no reason why the treatment should 
not apply to non-symmetric strains equally well. Note 
that the dependence of the strain on the stress is only 
through the responses of both these modes of motion to 
local stresses. This is, in fact, the main difference be- 
tween this theory and elasticity-based descriptions that 
relate directly the strain to the stress. 
The (symmetrised) strain due to rotation at the centroid 
of auxeton g is given directly by the tensor [2^ , 

= Q^oii (6) 

where is its angle of rotation, which depends on the 
local stress. Eq. ^ is written so that it holds both in 
2D and in 3D. In 3D, this expression is symmetric under 
exchange of the indices i and j, but anti-symmetric under 
exchange of k and I. This is due to the anti-symmetric 
nature of the description of the axes of rotations kl. In 
2D, there is only one axis of rotation, perpendicular to 
the plane, and the indices kl are redundant, which re- 
duces Qijki to the tensor Qj^ of eq. ([5]). This expression 
has been derived first in [25| for granular media, where 
it gives rise to dilation. It comprises the only relevant 
contribution to the strain when the auxetons are rigid 
and, as such should also describe well the systems dis- 
cussed in [20|. For what follows, it is important to note 
the observation in [3 that Q is a measure of the 
local rotational (or chiral) deviation of the auxeton from 
a global zero average. This rotation is best quantified by 
the sign of Tr{Q}. 

When elements can also fold and unfold, their expan- 
sions depend on the local stress. Significantly, there is no 
reason to expect that auxeton expansions be isotropic; 
depending on the choice of shape and the local structure 
around them, auxetons may expand differently in differ- 
ent directions. The expansive strain rate can be related 
directly to the local stress via 

where the non-isotropic expansion can be modelled into 
and different auxeton shapes would be described by 



different such matrices. Limiting the description to sym- 
metric strains imposes some constraints on the local ex- 
pansion tensor E, making its properties similar to those 
of the conventional compliance matrix in linear elastic- 
ity. However, such similarity would not exist for non- 
symmetric strains. For example, for such strains, E need 
not be symmetric under exchange of i and j. It is im- 
portant to note that, whilst the strain may have non- 
symmetric components, for example to describe large- 
scale vorticity, the stress cannot if there are no external 
couples to balance the excessive torque. This is one of 
the reasons that the following theory cannot be mapped 
readily to elasticity theory, nor to Cosserat theory. 
The total strain can be written then as 

^!j = %ki<i + QU(^kd'j^) (8) 

This relation is reminiscent of the yield equations in gran- 
ular systems (25|, but for two important differences. One 
is that, in granular systems, the rotating elements (the 
grains) can also slide relative to neighbours, a mechanism 
that auxetons do not possess, which gives rise to an ad- 
ditional, plasticity-like term. The other difference is that 
auxetons can fold and unfold (the _E-dependent term), 
which rigid grains cannot. 

Relation ([S]) makes good sense on the auxeton level. 
However, to be of use to materials that contain many 
auxetons, it must be coarse-grained (homogenised) to the 
continuum. To this end, one must average it over small 
volumes, containing sufficiently many auxetons. The ex- 
pansive term on the right hand side of ([SJ gives no prob- 
lems - one can average E and a independently to obtain 
a continuum-scale contribution. This is no different than 
the practice in conventional elasticity and plasticity mod- 
els. 

In contrast, the rotational term requires a careful con- 
sideration. Coarse-graining over the rotation field of a 
region of volume V, (9) = {i/V)J2g(^^^ can be carried 
out by replacing the volume average by a surface sum 
(or integral, for large enough regions), using Stokes the- 
orem. This leads immediately to the observation that 
the contribution to such an average comes only from the 
boundary of the region. Hence, if the system does not ro- 
tate globally, then the rotation per auxeton decays fast 
as the averaging volume increases and the macroscopic 
rotation has a zero average. It turns out that the tensor 
Q possesses exactly the same property. Since this tensor 
measure the local chiral fluctuation of an element, an av- 
erage over a region decays to the global zero average at 
exactly the same rate as {6). 

On the face of it, these two observations may seem to 
imply that the rotational contribution to the strain van- 
ishes on large scales. This, however, is not the case. The 
reason is that both 9 and Q possess the same local anti- 
correlations: when one auxeton rotates in one direction, 
elements in contact with it are more likely to rotate in 
the opposite, rather than in the same, direction. Simi- 
larly, if the tensor measures the rotation of an auxeton 
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at a given direction, nearest-neighbours of g are more 
likely than not to have Q's whose trace has the oppo- 
site sign. This anti-correlation means that, while each of 
these terms averages to zero independently over increas- 
ing regions, their product (QijkiOki) adds constructively 
over nearest neighbours, leading to a finite large-scale av- 
erage. It is exactly this average that leads to measurable 
bulk strain due to rotations of rigid particles in granular 
systems (dilation). We therefore conclude that eq. (H]) 
has a well defined homogenised large-scale version 



Bij = EijklCTkl + Qijkl&kl (f) (9) 

The only remaining question is how to derive local con- 
tinuous expressions for the rotational term. This can be 
done, using the renormalisation approach taken in (26| . 
A word of caution: the existence of a macro-scale con- 
tinuous theory does not imply that the strain is auxetic. 
Relation (j9|) gives the correct dependence of the strain 
on the local fields, but whether the ratio between per- 
pendicular strains (Poisson's ratio) is negative or positive 
depends on the relative contribution of the two terms on 
the right hand side of this relation. 

The advantage of relation ([9]) is that it identifies the 
precise role that the local structure plays in the coupling 
to the strain and the stress. As such, it is an important 
ingredient in the field equations of auxeticity theory. To 
complete the theory for I A structures, one still needs the 
local rotational response to the local stress, 9ki{cr). This 
relation is still missing and work to derive it is ongoing. 
Thus, the full set of auxeticity field equations in d dimen- 
sions consist of: 

(i) d{d +l)/2 balance eqs. (P)-©; 

(ii) d{d- l)/2 strain eqs. ©; 

(iii) d{d— 1)/2 rotation - stress response relations, Oki{cr). 
As in any theory, constitutive information is required. 
For the theory described here, this comprises the consti- 
tutive tensors E and Q, which could be obtained either 
phenomenologically or modelled theoretically for specific 
structures. 

The solution for quasi-static deformation then pro- 
ceeds as follows. First, one solves for the stress field 
from eqs. (II])-©. From this solution one finds the local 
rotational field Oki{cr), using the local rotation - stress 
relation. Substitution of the rotational field, the consti- 
tutive fabric tensor Q and the expansion tensor E into 
eq. ^ one then derives the total local strain. 

Since the dependence of the strain on local rotation 
and expansion of elements is valid regardless of the stress 



state, then all this theory, but for the stress field equa- 
tions, applies to any auxetic material. In particular, 
it applies to EA materials, where the stress equations 
should then be replaced by those of elasticity theory. In 
other words only the closure relation (|3]) is replaced by 
Saint tenant compatibility conditions [2 7|. supplemented 
with phenomenological or modelled expression for the 
stress-strain relations. 

IV. Discussion and conclusions 

To conclude, we have described a theory for strains in 
auxetic materials. The theory's main contribution is the 
explicit relation between the local auxetic strain and the 
local rotation and expansion of auxetons - the elementary 
building blocks of auxetic materials. This is a refinement 
of the currently existing elasticity theory which lumps 
these two contributions together into a stress-strain rela- 
tion. The identification of these strain mechanisms makes 
it possible to eventually arrive such a relation, since the 
local magnitudes of auxeton rotations and expansions do 
depend on the local stress. However, the explicit decom- 
position to rotation and expansion give insight into the 
correct symmetryies and details of such a stress-strain 
relation. 

Furthermore, using elasticity theory for lA would lead 
to erroneous results, which originate from two sources. 
Firstly, the stress state cannot be derived from elasticity 
theory and is likely to exhibit non-uniform force-chain- 
like fields. Secondly, the rotational and expansion re- 
sponses to the stress are of completely different nature. 
For example, the averaging properties of Q and E are 
completely different - while the latter has a well-defined 
macroscopic homogenised value, the latter does not. This 
is despite both terms having homogenised large-scale con- 
tributions. 

A significant implication of the above is that all aux- 
etic, whether lA or EA, must follow the universal strain 
eq. (jni). However, the stress state, which determines the 
local rotation and expansion of auxetons, depends on the 
correct stress description and this may vary between dif- 
ferent families of materials - isostaticity theory for lA and 
elasticity theory for EA. This then leads to the intrigu- 
ing conclusion that the auxetic behaviour of lA and EA 
materials should be markedly different, with the former 
exhibiting more non-uniform local auxetic behaviour. 

It is emphasised that eq. (0) does not ensure auxeticity, 
but rather it describes correctly the strain as a function 
of the local rotational and expansive fields. Whether the 
material exhibits a bulk negative Poisson's ratio depends 
on the different contributions of the two terms in the 
strain relation. 
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